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ABSTRACT 

Context. Oxygen is the third most common element in the Universe. The measurement of oxygen lines in metal-poor unevolved stars, 
in particular near-UV OH lines, can provide invaluable information about the properties of the Early Galaxy. 

Aims. Near-UV OH lines constitute an important tool to derive oxygen abundances in metal-poor dwarf stars. Therefore, it is impor- 
tant to correctly model the line formation of OH lines, especially in metal-poor stars, where 3D hydrodynamical models commonly 
predict cooler temperatures than plane-parallel hydrostatic models in the upper photosphere. 

Methods. We have made use of a grid of 52 3D hydrodynamical model atmospheres for dwarf stars computed with the code 
CO s BOLD, extracted from the more extended CIFIST grid. The 52 models cover the effective temperature range 5000-6500 K, 
the surface gravity range 3.5^.5 and the metallicity range -3 < [Fe/H] < 0. 

Results. We determine 3D-LTE abundance corrections in all 52 3D models for several OH lines and Fe i lines of different excitation 
potentials. These 3D-LTE corrections are generally negative and reach values of roughly -1 dex (for the OH 3167 with excitation 
potential of approximately 1 eV) for the higher temperatures and surface gravities. 

Conclusions. We apply these 3D-LTE corrections to the individual O abundances derived from OH lines of a sample the metal-poor 
dwarf stars reported in Israelian et al.(1998, 2001) and Boesgaard et al.(1999) by interpolating the stellar parameters of the dwarfs 
in the grid of 3D-LTE corrections. The new 3D-LTE [O/Fe] ratio still keeps a similar trend as the 1D-LTE, i.e, increasing towards 
lower [Fe/H] values. We applied 1D-NLTE corrections to 3D Fei abundances and still see an increasing [O/Fe] ratio towards lower 
metallicites. However, the Galactic [O/Fe] ratio must be revisited once 3D-NLTE corrections become available for OH and Fe lines 
for a grid of 3D hydrodynamical model atmospheres. 

Key words, nuclear reactions, nucleosynthesis, abundances - Galaxy:halo - Galaxy:abundances - Galaxy:evolution - stars: 
Population II - stars: abundances 



. _ 1 ■ Introduction There have been numerous studies of the oxygen abundance 

.,_,.,, in halo stars. Despite considerable observational and theoreti- 

^S , The metal-poor stars of the Galactic halo provide the fossil ^ efforts> ^ trend of [Q/Fe] ^ yersus [Fe/H] ig m un _ 

03 ■ record of the early Galaxy s composition. Dwarf halo stars are dear The ana , ig of ^ forbidden lme Ql 63QQ A in iants 

particularly relevant, because their atmospheres are not sigmf- shows a lateau ^ [Q/Fe] Q 4 _ Q 5 {n fce metallich 
icantly altered by internal mixing and provide a unique tracer _ 25 < [Fe/H] < _ { (Barbuy 19g8) and [Q/Fe] „ 

to constrain Galactic chemical evolutionary models. Oxygen is 7 for _ 4 < [Fe/H] < _ 2 .5 (Cavrel et al. 2004). A simi- 

a key element in this scenario, because it is the most abundant Jar behaviour is seen for metal . poor subgiant sta rs in the range 

element in stars after H and He. It is produced in the interiors of _ 3 < [Fe/R] < _ L5 ([0/Fe] „ QA _ Q ^ Garcfa Perez et al 

massive stars by hydrostatic burning and its content is modified 2QQ6 Ro ^ y ^ | Ga rcia P6rez et all d2006h already noted that 

during the explosive nucleosynthesis .in type II supernovae (SNe) by lotdn ^ measurements from the [OI] line for dwalfs 

and hypernovae (HNe) and returned into the interstellar medium (Mssen et al 2002)> sub iants (Garcfa Perez et al 2006) and { _ 

On the other hand, iron is created by both type II and type I SN ants (Cayrel et al 2004)> ^ icture ch s and an increasin 

explosions. However type I SNe progenitors have longer life- trend [Q/Pe] towards lowgr metallicities clearl 
times, which is why the abundance ratio of [O/Fe i] can be used 
to constrain the chemical evolution of the Galaxy. 
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The near infrared (I R) triplet Oi 7771-5 A in metal- 
poor dwarfs and subgiants ( Abia & Rebolo 1989; Israelian et al. 
2001; Garcia Perez et al. 2006) points towards increasing [O/Fe] 
values with decreasing [Fe/H], although the O abundances de- 
rived from the near-IR triplet are typically ~ 0.4-0.7 higher than 
those derived from the forbidden O i line dFulbright & Johnson] 
2003). The OH A-X electronic lines in the near ultraviolet 
(UV) provide also higher [O/Fe] ratios towards lower [Fe/H] 
value s in dwarf star s (Israelian et al. 1998; Boesgaard et al. 
Il999t lisraelian et a ll |200lk Gonzalez Hernandez et al. 2008). 
However, Garcia P erez et al.l (2006) found a quasi-plateau of 
[O/Fe] for subgiant stars in the range -3. < [Fe/H] < —1.5, us- 
ing Fe abundances determined from Fe n lines. The [O/Fe] ratio 
shows a negative slope if one instead uses the Fe abundances es- 
timated from Fe i lines, although with lower [O/Fe] values than 
those determined for metal-poor dwarfs. It is advisable to use 
Fe i lines instead of Fe n lines to derive the [O/Fe] ratio because 
of its similar sensitivity to the surface gravity. 

These O-abundance indicators present different complica- 
tions. The near-IR Oi triplet is susceptible to non-local ther- 
modynamical equilibrium (NLTE) effects (Kiselman 2001, and 
references therein), with abundance corrections below 0.2 dex, 
and is quite sensitive to the adopted r e jf. The [OI] is not sensi- 
tive to departures from LTE, but it is essentially undetectable 
in dwarfs with [Fe/H] < -2. The near-UV OH lines are 
stron gly sensitive to the temperature structure and inhomogini- 
ties ( Asplund & G arcia Perea 1200 U iGonzal ez Hernande z et al.l 
2008). In addition, Fe i li nes suffer from severe NLTE effects in 
metal-poor stars (see e .g. lThevenin & Id iart 1999). We note that 
Shchukin a et al.l (12005) have performed NLTE computations for 
the metal-poor subgiant HD140283 ([Fe/H] ~ - 2.5) with a sin- 
gle sn apshot of a 3D hydrodynamical simulation (Asplun d et al.l 
1999), and found NLTE-LTE corrections of +0.9 and +0.4 for 
Fe i and Fe n lines, respectively. 

Finally, the oxygen abundance in the Sun is still a mat- 
ter of debate. We will adopt throughout this work the value 
of log[JV(0)/iV(H] = 8.76, which was d etermined through 3- 
dimensional (3D) hydrodynamical models (Caf fau et al.l l2008). 

We have used a subset o f the CIFIST grid o f 3D hydrody- 
namical model atmospheres dLudwig e t al. 2009) to investigate 
the 3D-LTE and 3D-NLTE [O/Fe] and [O/H] trends in metal- 
poor dwarf stars. 

2. 3D hydrodynamical simulations 

The 3 D hydrodynamical model atmospheres (see iLudwig et al.l 
2009, for further details) were computed with the code 
CO s BOLD dFrevtag et al.ll2002k IWedemever et alj|2004l) . Each 
model consists of a representative set of snapshots sampling the 
temporal evolution of the photospheric flow. In Table Q] we pro- 
vide a summary of the 3D models used in this paper. The evolu- 
tionary time scale and the spatial scale of the convective surface 
flows roughly scale inversely proportional to the surface gravity. 
Following this scaling behaviour, we tried to ensure that the du- 
ration of each model sequence corresponds to at least 3600 s of 
solar-equivalent time. Moreover, we also scaled the size of the 
computational domain leading to about the same total number 
of convective cells in each model. 

The spatial resolution was typically (« x x n y x n z ) - 
(140, 140, 150), with a constant grid spacing in the (horizontal) 
x- and y-directions and a variable spacing in z-direction. Our 
standard solar model - actually not part of the grid used here, 
but mentioned for reference - has a horizontal size of 5600 km. 
The total height is 2250 km, from z = -1380 km below optical 



depth unity to z — 870 km above, located in the low chromo- 
sphere. The computational time step was typically 0.1-0.2 s. 




Fig. 1. 3D temperature structure of the 3D CO s BOLD atmo- 
spheric model, compared to the average temperature profile, 
(3D), and the IDlhd model. Darker shades indicate more likely 
temperatures. The stellar parameters and metallicity of this 3D 
simulation are r eff [K]/log(g[cm S - 2 ])/[Fe/H] = 5850/4/ - 3. 
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Fig. 2. Same as Fig. Q] but for the 3D model 

r eff [K]/log(£[cm s- 2 ])/[Fe/H] = 6270/4/ - 3. 



The description of the radiative energy exchange is impor- 
tant for the resulting temperature structure of a model. Here we 
provide some details about the binning-scheme which we ap- 
plied when modelling the wavelength dependence of the radia- 
tive transfer. The wavelength dependence is represented by 5 
multi-group bins in solar metallicity models and 6 bins at sub- 
solar meta llicities following the procedure origi nally laid out b y 
Nordlund (198 2]) and subs e quentl y refined bvlLudwig| (11 992); 
Ludwig. Jorda nT& Steffer] (11994 : IVogler. Bruls. & Schussler 
(2004). For test purposes we calculated a few models w ith 12 
bins. D ifferent from the original implementation of Nordlund 
dl982l) . the bin-averaged opacities were calculated explicitely 
from the run of the monochromatic opacity within a particular 
bin. No scaling among the bins was assumed. The sorting into 
wavelength groups was done applying thresholds in logarith- 
mic Rosseland optical depth {+oo,0.0, —1.5, -3.0, -4.5, -oo}for 
the 5-bin, and [+oo,0.1,0.0, -1.0, -2.0, -3.0, -oojfor the 6-bin 
schemes. For the 12 bins we used as thresholds {+oo, 0.15, 0.0, 
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-0.75, -1.5, -2.25, -3.0, -3.75, -4.5, -oo}; in addition, each of 
the first three continuum-like bins were split into 2 bins accord- 
ing to wavelength at 550, 600, and 650 nm. In all but one bin 
a switching between Rosseland and Planck averages was per- 
formed at a band-averaged Rosseland optical depth of 0.35; in 
the bin gathering the largest line opacities the Rosseland mean 
opacity was used throughout. The decisions about the number 
of bins, and sorting thresholds are motivated by comparing ra- 
diative fluxes and heating rates obtained by the binned opaci- 
ties in compariso n to high wav elength resolution. In comparison 
to the models of Asplund & Garcia Perez (2001 ), wh o worked 
with Stein-Nordlund models (Stein & Nordlund 1998) and used 
4 opacity bins, our treatment of the wavelength dependence of 
the radiation field puts extra emphasis on the continuum-forming 
layers. 

2.1. 3D temperature structure 

In Figs. Q] and [2] we depict the temperature structure of two 
3D hydrodynamical simulations for two different effective tem- 
peratures, 5850 and 6270 K, and the same surface gravity, 
log(g/cm s 2 ) = 4, and metallicity, [Fe/H] = -3. The main dif- 
ference between these two 3D models is the temperature inho- 
mogeneities which are larger in the hotter model. These strong 
temperature fluctuations are typically present for effective tem- 
peratures hotter than 5900 K and surface gravities lower than 4. 



Table 1. Details of the 3D hydrodynamical C0 5 BOLD model 
atmospheres. 




Fig. 3. Contribution functions to the equivalent width of the OH 
3 123 A line computed with C0 5 BOLD3D model (solid line) and 
IDlhd model (dashed line). The same is displayed for the OH 
3167 A line (3D model, dotted line, and IDlhd model, dashed- 
dotted line). The stellar parameters and metallicity of this simu- 
lation are T eS [K]/ log(g[cm S - 2 ])/[Fe/H] = 6270/4/ - 3. 

The 3D models displayed in Figs.Q]and[2]show an overcool- 
ing effect in the outer layers of the star with respect to ID mod- 
els. This is particularly relevant in all models with metallicities 
[Fe/H] = -3 and -2. 3D models with metallicities [Fe/H] = — 1 
and do not show a pronounced overcooling effect and the line 
formation is dominated by the temperature fluctuations. We note 
that the overcooling of the higher photospheric layers is the dom- 
inant 3D effect, but in the deep photosphere, at logr 1, 3D 

models are always slightly hotter than a corresponding hydro- 
static ID model in radiative-convective equilibrium. 

The comparison of 3D versus ID models depends on which 
particular ID model is chosen. We compared each of our full 
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fl Temporal average of the emergent T e g. 

b Total time covered by the temporal evolution of the photospheric 
flow named as snapshots. 



3D models to a corresponding (in T e g, logg, and metallicity) 
standard hydrostatic ID model atmosphere (hereafter denoted as 
IDlhd) and a ID model obtained from the temporal and horizon- 
tal average of the 3D structure over surfaces of equal (Rosseland) 
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Fig. 4. Same as Fig. |3]but for the Fe i 3849 A and 3852 A lines. 
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Fig. 5. Line profiles of the Fei 3849 A and 3852 A lines com- 
puted with C0 5 BOLD3D model (solid line for Fe i 3849 A 
and dotted line for Fei 3852 A) and IDlhd model (dashed 
line for Fei 3849 A and dashed-dotted line for Fei 3852 A). 
The stellar parameters and metallicity of this simulation are 
r eff [K]/log(g[cm s- 2 ])/[Fe/H] = 6270/4/ - 3. 



code Linfor3D0. It is necessary to assume a value of the micro- 
turbulence for the spectral synthesis of the ID models. We have 
adopted a £ t = 1 km s -1 . 

2.2. Line formation 

In Fig. [3] we display the contribution functions (CFs) 
to the equivalent width (EW) of the disc centre of the 
OH lines at 3123 A and 3167 A according to the 
C0 5 BOLD3D model with given stellar parameters and metal- 
licity r eff [K]/log(g[cm s- 2 ])/[Fe/H] = 6270/4/ - 3. These two 
features have different excitation potentials (\ = 0.2 and 1.1 eV 
for the OH 3 123 A and 3 167 A lines, respectively) and their con- 
tribution functions look similar in shape but different in strength. 
These molecular lines tend to form in the outer layers of the 3D 

model atmosphere, at log r between ~ -2.8 and 2.2, where 

the 3D model is cooler than the IDlhd model. Thus this line is 
significantly weaker in the IDlhd model and forms in the inner 
layers of the star. Because of the sensitivity of OH lines to the 
temperature, the line appears much stronger in 3D than in ID. 
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Fig. 6. 3D - IDlhd corrections for different T e ff and metallici- 
ties, computed for the OH 3167 A line. These 3D hydrodynami- 
cal model atmospheres have log g = 4. 



optical depth (hereafter denoted as (3D)). The IDlhd model 
is calculated with a ID atmosphere code called LHD. It as- 
sumes plane-parallel geometry and employs the same micro- 
physics (equation-of-state, opacities) as C0 5 BOLD. Convection 
is described by mixing-length theory. See Caf fau et alj £2007) 
for further details. Note that the IDlhd model depends on the 
choice of the mixing-length parameter, the (3D) model on de- 
tails of the averaging procedure. The temperature of the (3D) 
model is indeed obtained by averaging the fourth power of the 
3D temperature field. The choice was motivated by the property 
that this kind of averaging largely preserves the radiative flux 



that t his Kind ot ave raging largely 
(Steff en. Ludwig. & Frevtagffl 995). 



We adopted a mixing-length parameter a = 1 in the LHD 
models. However, this parameter only affects lines that are 
formed at optical depth log r > -0.5. As we will see in the next 
section, the OH lines in metal -poor ID models form at optical 
depth log t < - 1 . 

We compute the curves of growth of the OH lines and Fe i 
lines, for both ID and 3D models with the spectral synthesis 



This effect is u sually found in metal-poor atmospheres (see 
Asplund et al. 1999) and is particularly important for the oxygen 
abundances derived from OH molecules, as the differences are 
large in the layers where these lines are formed. The result is 
that the oxygen abundances become lower in the 3D formulation 
than in the ID. We quantify this effect through the 3D correction 
(see Sect. |231). In hotter 3D models (r eff > 5900 K) this effect 
becomes more severe, making the 3D corrections larger. 

We display the CFs to the EW of the Fe i lines at 3849 A 
and 3852 A in Fig. [4] These lines have different excitation po- 
tentials <x = 1.0 and 2.2 eV for the Fe 1 3849 A and 3852 A lines) 
and their CFs are different in both strength and shape. The Fei 

3849 A line extends from logr 4.5 to almost 0, whereas the 

Fe 1 3852 A line, due to the higher excitation potential, goes from 

logr 3.5 to 0. Nevertheless, the main contribution is located 

at logr ~ -3.5 for the Fei 3849 A line and at logr ~ — 1 for 
the Fe i 3852 A line. This means that the Fe i 3849 line is mostly 
formed in the outer layers of the 3D atmospheric model, and the 



2 more information on Linfor3D can be found in the following link: 
http://www.aip.de/~mst/Linfor3D/linfor_3D_manual.pdf 
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equivalent width is larger than in the IDlhd model. Although 
the Fe i 3852 A line is mainly formed in the inner layers, the 
EW computed for this 3D model is still larger than in ID, due to 
the contribution of the outer layers. Therefore, this line still pro- 
vides negative 3D corrections. In addition, although they have 
very similar excitation potentials, the ratio EW(3D)/EW( IDlhd) 



of the Fe 1 3849 A line is smaller than that of the OH 3 167 A line. 
This is probably caused by the stronger sensitivity of the OH line 
to the temperature structure of the 3D model. 

In Fig. [5] we display the line profiles of the Fei 3849 A 
and 3852 A lines for the same 3D and IDlhd models with 
stellar parameters and metallicity r e ff / log(g/cm s 2 )/[Fe/H] = 
6270/4/ - 3. The 3D profiles of both lines are more broadened 
than the IDlhd profiles in part because these ID profiles were 
not been convolved with a macroturbulent velocity, whereas in 
the 3D profiles the broadening effect due to macroturbulence is 
naturally included. In addition, the 3D profiles are slightly asym- 
metric due to the convective motions in the hydrodynamical sim- 
ulations included in the 3D model atmosphere. It is quite clear 
that the 3D profile of the Fe i 3849 A line is stronger than the 
IDlhd profile. This is not as evident for the Fei 3852 A line, 
where the line core is deeper in the IDlhd profile but the line 
wings are tighter. The ratios EW(3D)/EW( IDlhd) of these lines 
in these 3D and IDlhd models are 1.59 and 1.25 for the Fei 
3849 A and 3852 A lines, respectively. As a comparison, the ra- 
tios EW(3D)/EW( IDlhd) of the OH 3123 A and 3167 A lines 
are 20.2 and 8.2, respectively. Therefore, it is expected that all 
these lines should provide different 3D corrections. We will dis- 
cuss quantitatively these abundance corrections in Sect. 12.31 

2.3. 3D corrections 

The spectral synthesis code Linfor3D produces three curves of 
growth: full 3D, (3D) and 1D L hd dCaffau et all 120091 [2010). 
This allows us to estimate abundance corrections through the 
EW(3D) of the 3D model abundance to derive the correspond- 
ing abundances in the (3D) and IDlhd models. Two main effects 
distinguish 3D from ID models: the average temperature pro- 
file and the horizontal temperature fluctuations. We quantify the 
contribution of these two main effects by introducing the 3D cor- 
rection as 3D - IDlhd- The first and the second effect can also 
be separately estimated with the (3D) - IDlhd and 3D - (3D) 
corrections, respectively. The velocity fluctuations (i.e. intrinsic 
velocity fluctuations and microturbulence) are other effects that 
also play a role, but their influence on the derived abundance, al- 
though not negligible, is not as relevant. In Fig.|6]we display the 
3D - IDlhd abundance corrections of the OH 3167 A line for 
different effective temperatures and metallicities, and log g — 4. 
The 3D models with metallicities [Fe/H] = -1 and show very 
small and mainly positive 3D - IDlhd corrections, whereas 3D 
models with [Fe/H] = -3 and -2 show strong, negative correc- 
tions. The 3D - IDlhd generally increases towards hotter 3D 
models and lower metallicites. We note here that the behaviour 
of other OH lines is practically the same although for slightly 
different values for the 3D - IDlhd abundance corrections. 

In Figs.|7]and[8]we show the 3D - IDlhd abundance correc- 
tions of the Fei 3849 A and 3852 A lines for different effective 
temperatures and metallicities, and log g - 4. The behaviour of 
the Fe 1 3849 A line is quite similar to that of the OH 3 167 A line, 
although with smaller 3D - IDlhd values. On the other hand, the 
Fei 3852 A line behaves in a completely different way, and the 
3D - IDlhd corrections depend on whether the line formation 
occurs mainly in the outer or inner region of the atmosphere of 
a given 3D model at a given temperature and metallicity. 

In Fig. [9] we depict the (3D) - IDlhd corrections of the OH 
3167 A line for different effective temperatures and metallici- 
ties, and logg = 4. The (3D) - IDlhd corrections reflect how 
important these cooling effects are in the 3D models. Thus, the 
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Table 2. 3D abundance corrections for several OH lines using 3D CO s BOLD models with 6 and 12 opacity bins with the parameters 
r eff [K]/ log(g[cm s~ 2 ])/[Fe/H] = 6270/4/ - 3 and 6240/4/ - 3, respectively. 



Wavelength 
[A] 


X 

[eV] 


3D- 


- (3D),6-bin 

[dex] 


3D- 


- (3D).l2-bin 

[dex] 


A6-12bin 

[dex] 


3D- 


- lDLHD,6-bin 

[dex] 


3D- 


- lDLHD,12-bin 

[dex] 


A6-12bin 

[dex] 


3128.287 
3139.170 
3167.169 


0.210 
0.763 
1.108 




-0.852 
-0.647 
-0.544 




-0.953 
-0.735 
-0.622 


0.101 
0.088 
0.078 






-1.416 
-1.083 
-0.910 






-1.057 
-0.806 
-0.675 


-0.359 
-0.277 
-0.235 



strongest correction is found for the 3D model with metallicity 
[Fe/H] = -3 and r e ff ~ 5900. One also realizes by comparing 
Figs.|6]and|9]that the 3D models with metallicities [Fe/H] = -3 
and -2 and T e ff > 5900 have large temperature fluctuations that 
account for a significant fraction of the total 3D - IDlhd correc- 
tion. 

The sensitivity of the 3D - IDlhd corrections to the surface 
gravity is less relevant (see Fig. IT0T>. The derived 3D - IDlhd 
corrections are typically larger in 3D models with larger sur- 
face gravities, although only variations of < 0.2 dex are expected 
from logg = 3.5 to 4.5. 
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Fig. 11. 3D - IDlhd corrections for different excitation poten- 
tials and metallicities, computed for OH lines. These 3D hydro- 
dynamical models have T^ ~ 5900 K and log g - 4. 
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Fig. 10. 3D - IDlhd corrections for different surface gravities 
and metallicities, computed for the OH 3167 A line. These 3D 
hydrodynamical model atmospheres have T e g ~ 5900 K. 



The excitation potential determines the layer of the 3D 
model atmosphere where the line forms. Thus, it is expected 
that a line of the given molecule or atom forms further out if 
the excitation potential is lower. In the case of OH molecular 
lines, further out means lower temperature and therefore en- 
hanced molecule formation. This is exactly what it is displayed 
in Fig.[TT] the lower the excitation potential of a OH line is, the 
larger negative 3D - IDlhd corrections it provides. 

All in all, the dependence of the abundance corrections on 
metallicity shows a kind of threshold behaviour. Noticeable cor- 
rections indicating sizable temperature differences between ID 
and 3D mode l s set i n rather abruptly at metallicities below -1. 
iLudwig et al.l (2008) have argued that this can be expected on 
general grounds emerging from the temperature and metallicity 
sensitivity of the radiative relaxation time scale. 



2.4. Dependence on carbon abundance 

The molecular equilibrium of OH lines also involves the 
molecules CO, CH and C2. Thus, the C abundance surely af- 
fects the contribution function to the EW of the OH lines. 
The CO molecule is the most tightly bound and has the ten- 
denc y to get hold of all available O atoms, to the detriment of 
OH. iBehara et al.l (12010) have already checked this effect in C- 
enhanced metal-poor stars, showing that if the C abundance is 
large enough, the 3D - IDlhd abundance corrections become 
much smaller, reaching values of the order of -0.2 if the carbon- 
to-oxygen ratio is as high as [C/O] ~ 1 . In Table [3] we provide 
the 3D - IDlhd corrections to the O abundance derived from an 
OH line with excitation potential -0.1 eV, for several values of 
[C/O], which have been computed only for testing purposes with 
the 3D model with parameters T s g [K]/ log(g[cm s~ 2 ])/[Fe/H] = 
6270/4/ - 3. For these low ratios [C/O], the 3D - 1D LH d cor- 
rections are similar and within 0.15 dex. In addition, looking at 
these corrections in detail one realizes that (i) when [C/O] is low 
enough, the 3D - IDlhd correction to the O abundance derived 
from OH lines is only slightly dependent on the C abundance 
with changes < 0.05 dex, (ii) when the [C/O] is near the value 
zero, the C abundance starts to become relevant, with changes 
~ 0.20 dex, (iii) only when [C/O] is large enough, the 3D - 
IDlhd correction becomes negligible. 

In our computations we assumed that the carbon abundance 
scales with metallicity, [C/Fe] = 0. In models with metallicities 
lower than -1 dex, we have assumed the a-elements enhanced 
by a factor of [or/Fe] = 0.4. Therefore, the [C/O] = -0.4, which 
means that at the lowest metallicities, i.e. at [Fe/H] = -2 and -3, 
where the observed abundances can reach values of TO/F el ~ 
0.5 - 1.0 (see Sect. E) and [C/Fe] ~ (e.g. lLai et al.l l2007). our 
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estimated 3D - IDlhd corrections might have to be corrected 
downwards by 0.06 - 0.12 dex. This small effect has not been 
taken into account and goes in the direction of slightly increasing 
the strength of the 3D - 1D LH d corrections. 



Table 3. 3D - IDlhd O abundance corrections for several 
[C/O] ratios from the OH line 3106 A with an excitation po- 
tential equal to 0. 1 eV, using the 3D model with the parameters 
r eff [K]/ log(g[cm s- 2 ])/[Fe/H] = 6270/4/ - 3. 



Table 4. Selected lines for the abundance analysis. 



[C/O] 
[dex] 



[C/Fe] 
[dex] 



[O/Fe] 
[dex] 



3D 



IDlhd 
[dex] 



-1.0 
-0.6 
-0.2 



0.0 
0.4 
0.8 



1.0 
1.0 
1.0 



-1.53 
-1.49 
-1.37 



2.5. Dependence on opacity binning 

We tested a few 3D hydrodynamical models with 12 opacity bins 
and found that the cooling effect at log t < -2 is less pronounced 
than in our standard 6 bin models. From this result we expect 
that in the next generation of 3D models, using a refined opac- 
ity binning, the 3D - IDlhd abundance corrections would come 
out smaller for OH lines and Fe lines with excitation potentials 
below 1 .5 e V. Therefore, the corrections presented here for these 
lines can be considered as upper limits. 

In Table|2]we provide a comparison between the 3D - IDlhd 
corrections in the 6-bin and 12-bin 3D models. The 6-bin model 
produces larger corrections for these three OH lines, between ~ 
-0.36 and -0.23 dex, with excitation potentials between 0.2 and 
1.1 eV. In addition, the 12-bin model has stronger temperature 
fluctuations than the 6-bin model, since the 3D - (3D) is larger 
for the 12-bin model. This difference explains some of the larger 
scatter in oxygen abundances obtained from different lines when 
using 3D models than in the ID case (see Sect. [3]). 



3. Applications to the Galactic oxygen abundance 

The aim of this study is to apply our grid of the 3D - IDlhd 
corrections to a sample of metal-poor dwarf stars with avail- 
able O abundances derived from near- UV OH lines. This sam - 
ple was obtained from the studies of Boesgaa rd et alJ 11999), 
llsraelian et alj (1998, 2001), and the most metal-poor dwar f stars 
of the binary CS 22876-032 from Gonzalez Hernandez et al. 
(2008). We adopted the stellar parameters and the O and Fei 
abundances of the dwarf stars provided in those works. In Fig.[T3l 
we display the 1D-LTE [O/H] and 1D-NLTE [O/Fei] ratios ver- 
sus 1D-NLTE [Fei/H] in metal-poor dwarf stars. 

We applied the same formula given in Israelian et al. (2001) 
to correct the ID Fei abundances for NLTE. This formula was 
derived by these authors by fiting the ID NLTE-LTE corrections 
provided by Thev enin & Idia rt (1999) as a function of metal- 
licity. They also noted that these NLTE corrections were not 
significantly sensitive to the stellar parameters in their sample 
of metal-poor stars. This formula provides corrections of +0.25 
and +0.37 at [Fe/H] ~ -2 and -4, respectively. In the liter- 
ature there are other studies of the NLTE effect on Fe lines 
(e.g. iKorn et al]|2003t IShchukina et aT 2005: M ashonkina et alj 
2010) in metal-poor stars, but none of them provides a similar 



Wavelength 


Transition 


A 


logg/ 


[A] 




[eV] 




3123.949 


OH/>„(9.5) 


0.205 


-1.960 


3128.287 


OH P 2 2(8.5) 


0.210 


-2.021 


3139.170 


OH 222(17.5) 


0.763 


-1.559 


3167.169 


OH 2^22(21.5) 


1.108 


-1.541 


3255.493 


OH P 22 (23.5) 


1.299 


-1.809 


3843.257 


Fei 


3.047 


-0.241 


3849.967 


Fei 


1.011 


-0.871 


3852.573 


Fei 


2.176 


-1.185 



table with NLTE-LTE corrections for a grid of models at dif- 
ferent meta llicities. NLT E effects on the OH lines may be also 
significant. lAsplund & Garcia Perez (2001) tried to investigate 
this effect using a two-level OH molecule and found NLTE cor- 
rections up to +0.25 dex at metallicities below -2. We note here 
that in this study we have not taken into account any possible 
NLTE effect on the OH lines. The trend displayed in Fig. [13] 
clearly shows an increasing [O/Fe i] ratio towards lower metal- 
licities. The error bars in this figure only show the dispersion of 
the abundance measurements using ID models. We derived the 
slope of a linear fit to all data points shown in the lower panel 
in Fig. [13] providing -0.39 ± 0.03. This value ma y be compared 
the slo pe of the 1D-NLTE [O/Fe] ratio reported in Israelian et alJ 
(2001), whose value is -0.33 + 0.02. The d ifference is likely 
related to the fact that Israelian et al. (2001) used data points 
from Israelian et al. (1998, 2001), Boesgaard et al. (1999) and 
lEdvardsson et al.l (ll993l) . 

We determined for all 3D models in Table Q] the 3D - 
IDlhd abundance corrections for the OH and Fe i lines given in 
Table [4] The Fe i lines have three different excitation potentials 
from roughly 1 to 3 eV. W e chose some stars in the sample of 
llsraelian et al.l d 1 998l 1200 lh . for which we have near-UV spectra 
to check if the selected Fe i lines provide appropriate [Fe/ H] val- 
ues. For this estimate we used the code ID LTE MOOG (Sneden 
1973) and LTE model atmospheres with q-ele ments enhanced, 
by +0 .4 dex, computed w ith the Linux version ( Sbor done et al.l 
120041) of the ATLAS code dKuruczll993h. We used the new opac- 
ity distribution functions (ODFs^ of ICastelli & Kuruczl (l2003h 
with the corresponding metallicity. In Table [5] we give the ID 
LTE Fe abundances from these Fe i lines. The Fe i 3849 A line is 
quite strong and starts to be in the saturation regime for metal- 
licities above -2. If one takes into account that the 3D-LTE cor- 
rections for Fe i lines of different excitation potentials are dif- 
ferent, one should obtain different abundances when using ID 
models, according to Fig.QT] We hardly see this strong effect in 
the values given in Table[5] In Fig. [12] we depict these Fe i abun- 
dances for different values of microturbulent velocities for three 
of the stars in Table[5] As soon as one goes to higher metallicities 
and lower excitation potentials, the Fe i is more sensitive to the 
adopted microturbulence. However, at the lowest metallicities, 
the Fe i line at 3 eV provides more uncertain results because the 
EW is very low and therefore more sensitive to the quality of the 
observed spectra. We believe that our choice of microturbulent 
velocity equal to 1 km s is a reasonable assumption, although 
the hottest stars, with T e s ~ 6300 K, in our sample may require 
a slightly higher microturbulence. 

In order to derive 3D abundances, we interpolated within the 
grid of 3D - IDlhd corrections using the stellar parameters and 
metallicities of the sample stars. For OH lines, we applied indi- 
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Fig. 12. ID LTE [Fe/H] abundances of some the stars in Table |5] Fig. 14. 3D-LTE [O/H] and 3D-NLTE [O/Fe i] ratios versus 3D- 
excitation potential of the Fe i line, for three different microtur- NLTE Fei abundances computed with 3D-LTE models. The 
bulent velocities adopted in each computation. symbols are the same as Fig. [13] 



vidual abundance 3D - IDlhd corrections to the O abundance 
derived from each OH line and finally computed the average 3D 
O abundance of all available lines in each star. The error of the 
[O/H] abundance is then estimated from the dispersion of the 
3D-LTE abundances of all OH lines. 




-3 -2 -1 

[FeI/H],„ K . L , b 

Fig. 13. ID [O/Fei] ratio versus ID Fei abundances in 
metal -poor dwarf stars computed with 1D-LTE models us- 
ing OH lines from Boesgaard et al.(1999, triangles), Israelian 
et al.(1998, 2001, diamonds). The circles represent the most 
metal-poor dwarf stars of the binary CS 22876-032 from 
iGonzalez Hernandez etal] d2008l) . The stars labeled "A" and "B" 
represent the primary and secondary star in this binary. 



For Fei we us ed a different recip e . The papers by 
llsraelian et"ai] d!998h and iBoesgaard et al] d 19991) do not pro- 
vide the Fe i lines used to determine the metallicity. On the other 
hand, the paper by Israe lian et al.l (2001 ) gives a few Fe i lines in 
the near-UV spectral region where the OH lines form. The exci- 
tation potential of these Fe i lines is^ ~ 1 eV. However, these Fe i 
line are too strong and gets saturated in stars with [Fe/H] > -2. 
Therefore, we decided to adopt the following strategy: (i) for 
stars with [Fe/H] < -2 we only used the Fei 3849 A line to ap- 
ply a 3D - IDlhd correction to the ID metallicities; and (ii) for 
stars with [Fe/H] > -2 we derived an average correction with the 



two other Fe i lines at 3843 and 3852 A. For [Fe/H] > -2, we es- 
timated the error of the 3D-LTE [Fe i/H] by adding quadratically 
the dispersion of the 3D - IDlhd corrections and the dispersion 
of the 1D-LTE abundances. For [Fe/H] < -2, we adopted ad hoc 
a dispersion of 0.3 dex in the 3D - IDlhd abundance correc- 
tions because we were only using one Fei line. In Sect. 12.31 we 
discussed the expected differences in the 3D-LTE 3D - IDlhd 
correction between the Fe i 3849 A and 3852 A lines for models 
with [Fe/H] < -2. Thus, it is clear that if one would use both 
Fei lines to estimate the 3D - IDlhd correction for stars with 
[Fe/H] < -2, one would get smaller corrections, but the error 
bar associated to these corrections would be very large. 

As a first approximation, we corrected our values for the 3D- 
LTE Fei abundances for NLTE effects with the 1D-NLTE cor- 
rections already applied before to the ID case. Finally, we deter- 
mined the 3D-NLTE [O/Fe i] ratio by computing the difference 
between the 3D-LTE O abundance and the 3D-NLTE Fe i abun- 
dance in all dwarf stars of the sample. We are aware that these 
NLTE corrections may be too small, but we also know that our 
3D-LTE corrections in both OH and Fe lines may be also overes- 
timated because of the opacity binning used to compute our grid 
of 3D models (see Sect. 12.5b . In addition, the carbon abundance 
may be another factor to take into account (see Sect. l2.4l i. We see 
that these effects tend to go in opposite directions when trying to 
define the slope, if it exists, of the ratio [O/Fe] versus [Fe/H] at 
the lowest metallicities. 

In Fig. d we depict the 3D-LTE [O/H] and 3D-NLTE 
[O/Fe i] trend in metal-poor dwarf stars. The slope of the linear 
fit to the 3D-NLTE [O/Fe i] trend is smaller in absolute value, 
-0.16 + 0.04, but still negative and relatively significant, at least 
under the assumptions made in this work. The error bars are 
larger in this figure than in Fig [13] because we added quadrat- 
ically the standard deviation on the OH and Fe abundances, and 
the dispersion on the 3D - IDlhd abundance corrections from 
different OH and Fe lines. 

There is also a larger scatter in the almost linear relation be- 
tween the 3D-NLTE [O/Fe] and [Fe/H] ratios. The reason may 
be that the 3D corrections applied to the stars depend not only 
on the metallicity of the star, but also on the stellar parameters. 
Therefore, when we apply the 3D - IDlhd corrections, differ- 
ent stars move to different directions in the plane [O/Fe] versus 
[Fe/H]. As we stated in Sect. \T\ large and positive 3D-NLTE 
corrections for Fei lines are expected in metal-poor 3D mod- 
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Table 5. ID Fe abundances for six stars in the sample of llsraelian et al] (Il998ll200lh . from three Fei lines of different excitation 
potentials 



Star Ta \ogg £ [Fe/H] [Fei/H] 3 i 43 [Fei/H] 3 i 4 9 [Fei/H] 3152 

[K] [dex] [km s~' ] [dex] [dex] [dex] [dex] 



G 64-12 


6318 


4.20 


1.0 


-3.37 


-3.44 


-3.38 


-3.43 


G 64-37 


6310 


4.20 


1.0 


-3.22 


-3.22 


-3.21 


-3.19 


LP 815-43 


6265 


4.54 


1.0 


-3.05 


-3.06 


-2.92 


-3.05 


HD 140283 


5550 


3.35 


1.0 


-2.70 


-2.77 


-2.57 


-2.83 


HD 166913 


5965 


3.90 


1.0 


-1.80 


- 


-1.57 


-1.78 


HD 76932 


5800 


3.85 


1.0 


-1.10 


- 


-1.48 


-1.45 



els dShchukin a et al] 120051) . but this needs to be done for the 
whole grid of 3D hydrodynamical models to see how the pic- 
ture changes. We also note here that the computations done by 
Shchukinaet al. (2005) do not consider the inelastic collisions 
with neutral hydrogen atoms, because there is at present no re- 
liable values for the collisional rates, and the classical Drawin 
formula leads to uncertain estimates. They also pointed out that 
inelastic collisions with neutral hydrogen atoms would tend to 
alleviate the NLTE effects, and therefore, the difference between 
the NLTE and LTE may be considered as the maximum ef- 
fect that inelastic collisions might produce. Figure Q3] shows the 
3D-NLTE [O/Fe i] ratio versus effective temperature of the star. 
There is no clear trend with effective temperature, except that 
due to selection effects, most of the dwarf stars with higher 3D- 
NLTE [O/Fe i] values are hot stars (T e s > 5900), with the excep- 
tion of very few cases. 

Finally, we remark here that after all these corrections, the 
trend in [O/Fe i] obtained from near-UV OH lines in metal-poor 
dwarf stars is consistent with [O/Fe] ratios obtained from the for- 
bidden [O i] in metal-poor giant star s for which the 3D correc- 
tions are expected to be negligible (Gonzalez Hernandez et al. 
2008). 



Israelian et al. (1998,2001) 
Boesgaard et al. (1999) 




Fig. 15. 3D [O/Fe i] ratio versus effective temperature. The sym- 
bols are the same as Fig. [13] 



4. Summary 

The large grid of 3D hydrodynamical model atmospheres of 
dwarf stars quite captivate d our attention during the last three 
years dLudwig et al.ll2.Q09h . It took thousands of hours in com- 
puting time to build such a grid. We have used 52 3D models 



extracted from this grid with the following stellar parameters 
and metallicities: 7/ eff = 5000, 5500, 5900, 6300 and 6500 K, 
logg = 3.5, 4 and 4.5, and [Fe/H] = 0, -1, -2 and -3. The 
main difference with the "classical" ID models is that 3D mod- 
els show temperature inhomogeneities and a cooler average tem- 
perature profile in models with [Fe/H] < — 1. 

This allowed us for the first time to compute 3D abundance 
corrections of several near-UV OH and Fe i lines. These 3D cor- 
rections are generally larger for higher effective temperatures, 
larger surface gravities, and lower metallicities. In addition, lines 
with lower excitation potentials show stronger 3D corrections. 

We applied this grid of 3 D corrections to a sampl e 
of metal- poor dwarf s tars from Isra elian et alJ (I1998. 2001), 
Boesgaar d et alj ([1999) and the most metal-po or dwarf stars 
of the binary CS 22876-032 from Gonzalez Hernandez et al. 
(2008). We interpolated within this grid, using the stellar param- 
eters and metallicities of these stars. 

Finally, we are able to display a new trend with the 3D-NLTE 
[O/Fe i] ratio versus metallicity and this trend still increases to- 
wards lower metallicities, as the 1D-NLTE, but with a smaller 
slope in absolute value. However, we caution that some assump- 
tions (as e.g. the restriction to a 6-bin scheme in most 3D models 
and NLTE corrections only in ID), which we believe to be rea- 
sonable, have been adopted to achieve this result and that there- 
fore, this result should be taken with that in mind. A full 3D- 
NLTE study of all 3D model atmospheres presented in this paper 
must be performed for both near-UV OH lines and Fe i lines to 
see if this trend changes. 
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